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Abstract 

A  family  of  integer  programs  is  considered  whose  right-hand-sides  lie 
on  a  given  line  segment  L.   This  family  is  called  a  parametric  integer 
program(PIP) .   Solving  a  (PIP)  means  finding  an  optimal  solution  for  every 
program  in  the  family.   It  is  shown  how  a  simple  generalization  of  the 
conventional  branch-and-bound  approach  to  integer  programming  makes 
it  possible  to  solve  such  a  (PIP) .  The  usual  bounding  test  is  extended 
from  a  comparison  of  two  point  values  to  a  comparison  of  two  functions 
defined  on  the  line  segment  L.  The  method  is  illustrated  on  a  small  example 
and  computational  results  for  some  larger  problems  are  reported. 
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1.  Introduction 

The  purpose  of  this  paper  Is  to  show  how  a  simple  generalization  of 
the  conventional  branch-and-bound  approach  to  integer  programming  makes 
it  possible  to  solve  a  parametric  Integer  program.   Following  Nauss  [6] 
we  shall  call  the  family  of  programs 


e' 


(P^)      max     I     r  X 
j=l  ^   -' 


subject   to 
n 


l^Vi  '  ^^''^i 


t 


"i 


Isl^m 


X     e   {0,1} 


ISj^n 


for  0<e<l  a  single  parametric  integer  program  (PIP) .      By  "solving"   (PIP) 

we   shall  mean  obtaining  an  optimal  solution   of    (P.)    for  every  0^9Sl   for  which 

(P  )    is   feasible.     We  assume  that    (P.)    is   feasible  for  at   least  one  value  of   i 

e  o 

Parametric  Integer  programming  has  only  recently  emerged  as  a  topic 
of  research.   The  pioneering  papers  include  Noltemeier  [7],  Roodman  [9«10], 
Piper  and  Zoltners  [8],  and  Bowman  [1].  Nauss  [6]  has  reviewed  this  earlier 
work,  and  contributed  many  new  results  for  parameterizatlons  of  the  objective 
function.   The  present  paper,  which  has  grown  out  of  the  authors'  work  on 
synthesizing  dynamic  programming  with  branch-and-bound  [3,4,5],  is  devoted 
to  the  right-hand-side  case. 
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In  parametric  linear  programming,  the  first  step  Is  to  solve  (Pf^)  . 

I.e.  (P„)  for  6=0.   Then  the  direction  vector  d-(d, d  )  Is  specified 

o  1       m 

and  the  analysis  Is  performed  by  driving  6  from  0  to  1.   Critical  values  of 
6  and  new  optimal  solutions  are  identified  one  at  a  time  as  6  increases. 
In  the  procedure  for  parametric  Integer  programming  to  be  presented  here, 
the  direction  d  must  be  specified  in  advance.   The  (PIP)  is  solved  in  one 
branch-and-bound  search.   The  usual  bounding  test  is  modified  so  that  a 
partial  solution  is  eliminated  only  if  none  of  its  descendants  is  optimal 

for  any  (P„) ,  O^e^l.   This  means  that  some  partial  solutions  must  be 

o 

retained  that   could  otherwise  be  eliminated  if  only   (P^)   were  of  Interest. 
The  severity  of  the   resulting  computational  burden   depends   on   the  magnitude 
of  d. 

The  organization  of  the  paper  is  as   follows.      A  prototype  branch- 
and-bound  algorithm   for   (p   )    is  presented  In  Section  2. 

The  lower  bound  and  upper  bound  functions   are  developed  in  Sections   3  and 
4,   respectively.      The  modified  branch-and-bound  algorithm  for  (PIP)    Is 
given  in  Section  5   and  applied  to  a  sample  problem  in  Section  6.    Computational 
experience  with  the  algorithm  is  reported  in  Section  7. 


2,     A  prototype  branch-and-bound  algorithm 

We  shall  draw  upon  the   framework  and  terminology  of  Geoffrlon  and 
Marsten    [2]    to  describe   a  simple   linear  programming  based  branch-and-bound 
algorithm  for   (Pfv)  •     Problem  (?_)    is  separated,  by   fixing  variables  at   zero 


I 
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and  one,  into  smaller  candidate  problems  (CP^) .  Each  candidate  problem 
has  an  associated  set  of  fixed  variables  F  £  J"{1,  ...,  n}  and  partial 
solution  X  .  That  is,  (CP*^)  is  defined  by  the  conditions  x  -  x^  for  j  e  F^. 
The  current  set  of  candidate  problems  is  called  the  candidate  list.   If 
any  feasible  solutions  of  (?»)  are  known,  the  best  of  these  is  called  the 
incumbent  and  its  value  denoted  by  LB.   If  we  let  J  =  J-F^  be  the  set  of 
"free"  variables  and 

jeF^  -'  -^ 

where  A.  is  the  jth  column  of  A,  then  a  typical  candidate  problem  may  be 
written  as 

(CP^)   I  r  x^  +  max  y    r  x 


subject  to      y    a.  .X.  s;  b.-3; 
5,jq  ij  J    i  1 


isism 


X   €{0,1}  JEj** 

An  upper  bound  on  the  value  of  (CP  )  is  obtained  by  solving  its  LP 

relfixation  (CP^  ) .   It  is  also  helpful  to  compute  a  lower  bound  on  the  value  of 

(CP^) .  This  can  be  done  by  using  a  heuristic  to  find  a  feasible  solution 

of  (CP^) .   This  feasible  solution,  if  it  is  better  than  the  incumbent, 

becomes  the  new  incumbent,   a  prototype  branch-and-bound  algorithm  may  now 

be  described  as  follows. 
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Step  1.   Place  (P.)  in  the  candidate  list  and  set  LB  -  -  «  . 

Step  2.   If  the  candidate  list  is  empty,  stop.   If  there  is  an 

incumbent,  it  is  optimal  for  (Pq).   Otherwise  (P_)  is  infeaalble. 

Step  3.   Select  a  candidate  problem  and  remove  it  from  the  candidate 
list.   Call  it  (CP''). 

Step  A.   Solve  the  linear  program  (CP^  ).  Let  Ub'^  denote  its  optimal 
value. 

Step  5.   If  UB*^  5  LB,  go  to  Step  2. 

Step  6.   If  the  optimal  solution  of  (CP^  )  is  all  integer,  make  this 
solution  the  new  incumbent;  set  LB  =»  UB  ,  and  go  to  Step  21 


Step  7.   Use  a  heuristic  to  find  a  feasible  solution  of  (CP  ).   Let 
H^  denote  its  value.   If  h''  >  LB,  ' 
the  new  incumbent  and  set  LB  ■  H  . 


H  denote  its  value.   If  H^  >  LB,  then  make  this  solution 


Step  8.   Separate  (CP^)  into  two  new  candidate  problems  (CP'^  )  and  (CP^  ) 

by  choosing  peJ^  and  setting  f'^  =  f'^   -  F*^  u{p},  x*^  -  0,  x''   -l. 

Place  (CP^  )  and  (CP^  )  in  the  candidate  list  and  return  to 
Step  2. 


A  great  many  variations  on  this  pattern  are  described  in  [2],  but  this 

prototype  will  suffice  for  our  purposes.   Step  5  is  the  bounding  test.   If 

this  test  is  satisfied,  then  no  descendant  of  x'^  is  better  than  the  incumbent. 

Notice  that  the  bounding  test  includes  the  case  where  (CP  ) ,  and  hence 

K 

(CP*^)  ,  is  infeasible  since  then  UB^  =  -  ».   if  (CP^)  does  not  have  to  be 
separated  at  Step  8,  then  we  say  that  it  has  been  fathomed.   This  occurs 


if  (CP'^)  passes  the  bounding  test  or  if  (CP^  )  has  an  all  integer  solution. 
Step  7,  the  heuristic,  is  optional.  Its  purpose  is  to  strengthen  the  bounding 
test  by  inqsroving  the  incumbent  and  increasing  LB. 


•  '^- 


The  modifications  that  must  be  made  to  this  prototype  algorithm  to 
solve  (PIP)  are  confined  to  Steps  5,  6,  and  7.   The  notion  of  the  incumbent 
must  be  generalized  from  a  single  value  LB  to  a  function  LB (6)  defined  on 
0  ^6^1,   The  upper  bound  must  also  be  expressed  as  a  function  of  6:  Ub"(9), 
The  bounding  test  then  becomes  a  conparlson  of  two  functions  on  the 
interval  0^e<l  rather  than  just  a  point  comparison  for  6"0. 


3.   The  optimal  return  and  lower  bound  functions. 

In  this  section  we  shall  investigate  the  behavior  of  the  optimal 
value  of  an  Integer  program  as  a  function  of  its  right-hand-side.  Let 
the  optimal  return  function 

f(b')  =  max  rx    subject  to  Ax  s  b^ 
X  6  {0,1} 
be  defined  for  b'  e  R™.   It  is  apparent  that  f(b')  is  nondecreasing  in  each 


comp 


onent  of  b'.   Let  {x   I  k  e  K}  be  the  set  of  all  feasible  solutions  of 


(PIP),  i.e.  of  all  (P.)  for  0^651.   For  each  k  e  K,  define  the  step  function 

9 


f^(b 


')  -  j  j-i  ^  ^        j-i  ^  ^ 

L  _  oo   otherwise 


,m 


otherwise 
for  all  b'  e  R"*. 

The  optimal  return  function  f(b')  is  the  pointwise  maximum  of  this  finite 

collection  of  nondecreasing  step  functions 

f(b')  =  max  {f^(b')  |  k  e  K} 
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and  is  therefore  itself  a  nondecreasing  step  function. 

Now  suppose  that  the  Solutions  {x  |  k  e  K}  are  known,  where  K  £  K. 
A  lower  approximation  of  f(b')  may  be  constructed  from  these  known  solutions, 
namely 

f(b')  =  max  {f^(b')  I  k  €  K}. 

Clearly  fCb")  is  also  a  nondecreasing  step  function  and  is  a  lower  bound 
function  for  f(b'),  i.e.  f(b')  s  f(b')  for  all  b'eR"*.   The  approximation 
can  be  improved  as  new  feasible  solutions  are  discovered. 

We  are  interested  in  a  particular  "slice"  of  f(b')  and  f(b'):   the 
line  segment  {b+6d|0^e^l}  where  b  is  the  right-hand-side  of  (P.)  and  d  is 
the  given  direction  vector.   Define  g(e)-f(b+ed)  and  LB(e)=f (b+0d)  for 
0^e<l.   Then  g(e)  and  LB(e)  are  both  step  functions  and  LB(9)  5  g(e)  for  all 
Ose<l.   If  d>0,  then  g(e)  and  LB(e)  are  both  nondecreasing.  (See  Figure  1). 
There  is  at  least  one  optimal  solution  of  (PIP)  associated  with  each  step 
of  g(6).   Solving  (PIP)  is  equivalent  to  constructing  g(6)  by  finding  at 
least  one  x  solution  for  each  of  its  steps. 

The  procedure  for  constructing  LB(9)  from  the  known  feasible  solutions 
is  as  follows.   For  each  k  e  K  define 

,  n     , 

9^  -  min  {e  |  I     A  x^  s  b  +  9d}  (3.1) 

k  "     k 

9,  =  max  {9  I  j;  A  x^  s  b  +  9d}  (3.2) 

k    k 
where  9^  =  6  =  +  «>  if  the  indicated  set  is  empty.   Then 

LB^(e)=(L  ^.^^  if  9^<9^9^  (3.3) 

*■_  oo  otherwise 
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g(e) 

LB(e) 

i 


Figure  1.   Typical  g(e)  and  LB(e)  functions. 


Figure  2.   Typical  UB^(e)  and  UB*^(e;*)  functions. 


and 


k 
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LB(e)  -  max{LB  (6)  |  k  c  K}.  (3.4) 

The  solutions  which  determine  LB(e)  will  be  called  the  Incumbents.   Each  one 
is  Incumbent  for  a  particular  interval  of  8. 


4*  The  upper  bound  functions. 


Consider  a  given  partial  solution  x  .  In  order  to  demonstrate 
that  no  descendant  of  x  could  be  optimal  for  any  (P.) ,  we  need  an  upper  bound 
on  the  return  achieved  by  any  descendant  and  this  upper  bound  must  depend 
on  6.   Such  an  upper  bound  can  be  obtained  by  introducing  (ed)  into  the 
relaxed  candidate  problem  (CP^  ) .   Define 


UB^(e)  «  I         r,x^  +  max  y    r^x. 


subject  to 


iSlSm 


•.^...' 


so  that  Ub'^CO)  =  UB*^.   It  is  well  known  that  UB*^(6)  Is  concave  and  plecewlse 
linear  on  0<esl.  The  function  UB''(9)  could  be  obtained  explicitly  by 
ordinary  parametric  linear  programming.   The  computational  burden  Involved 
in  doing  this  for  every  candidate  problem  Could  be  quite  substantial,  however. 
Fortunately  any  dual  feasible  solution  of  (CP^  )  can  be  used  to  construct 
a  linear  upper  bound  fimction  for  UB^(e).   An  optimal  dual  solution  of 
(CP^  ),  barring  degeneracy,  yiiilds  the  first  linear  segment  of  UIJ^(O). 

K 
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By  linear  programming  duality  we  know  that; 

m  n 

m'^(e)   =  I    ^  T  x\  +  m±nl     u  (b  +ed  -bJ)  +  I      v 


subject   to 

i,  "i^j "  "j  ^  'j 

J.J" 

u^  s  0 

ISiSm 

v^^O 

l^j^n 

For  notational  convenience  we  have  included  all  of  the  v,  variables, 

even  though  v  =0  for  jeF^  in  any  optimal  solution.  Let  d"  denote  the 

dual  feasible  region 

m 
D^  =  {(u,v)^0  I  I       u  a   +  v^r  for  jeJ^}. 

Since  the  primal  variables  are  all  bounded  and  at  least  one  (P  )  is 
feasible,  we  may  conclude  that  D  is  non-empty.   Let  {(u  ,  v  ) | teT  }  and 
{(y^,  z^)|seS*^}  denote  the  sets  of  extreme  points  and  extreme  rays, 
respectively,  of  D^.  Taking  e=(l,  ...,  1)  we  have 

Ub'^CQ)  ^  I    r  x^  +  u^Cb+ed-e'^)  +  v'e 
jeF^  J  J 

for  all  tel'^,  with  equality  if  (u  ,  v  )  is  optimal  for  the  objective 

function  uCb+ed-B^^)  +  ve.   As  a  fimction  of  9  then,  the  return  achieved  by 

any  descendant  of  x  is  bounded  above  by: 

UB'^(e;t)  =  (u'^d)e  +  [  I         r  x\  +   u'^Cb-g'^)  +  v^e] 

jeF^  J  J 

for  any  teT^.   This  is  a  linear  function  of  6  and, since  u  sO,  it  is 

nondecreasing  if  d>0. 
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In  the  modified  branch-and-bound  edgorithm  for  (PIP) ,  linear 
programming  Is  applied  to  (CP^  )  as  usual.   The  functions  UB''(e;t)  are 
obtained  at  no  extra  cost.   The  function  obtained  from  an  optimal  dual 
solution  will  be  denoted  UB^(e;*).  Barring  degeneracy,  UB^(e;*)  coincides 
with  the  first  linear  segment  of  UB^(e).   (See  Figure  2).  As  in  conventional 
branch-and-bound,  if  the  dual  simplex  method  is  used,  then  suboptimal  dual 
solutions  can  be  used  to  perform  additional  weaker  tests. 


If  (CP^  )  is  infeasible,  then  the  simplex  method  will  terminate 

t     t  S    8 

with  an  extreme  point  (u  ,  v  )  s  0  and  an  extreme  ray  (y  ,  z  )  ^  0, 
such  that 

y^(b-B'^)  +  z^e  <  0. 

If  y^d^O,  then  we  may  conclude  that  UB^(e)  -  -  »  for  all  Osie^l.   If  y'^d>0, 
then  UB'^(e)  =»  -  «  for  O^eO*  and  UB^(e)  ^  UB^(e;t)  for  8*^6^1,  where 

e*  =  -y^Cb-e"^)  -  z^e 


y^d 


I 


5.     A  branch-and-bound  algorithm  for   (PIP). 


Now  that   the  upper  and  lower  bound   functions  have  been  derived, 
the  generalized  bounding  test  may  be   stated.      The  partial  solution  x'   doew 
not  have   a  descendant   that  is  better   than  an  incumbent    if 
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UB^(e)  <;  LB(9) 


or  if 


UB^(e;t)    <.  LB(e) 


for  all        O^esl 


for  all       O^e^l 


for  some  teT^.      This   test   is   the  basis   for  a  modified  branch-and-bound 
algorithm  that   can  solve   (PIP) . 


Step  1. 
Step  2. 

Step  3. 

Step  4. 

Step  5. 


Step  6. 
Step   7. 

Step   8. 


Place   (Pq)    in  the  candidate  list   and  set  LB(e)    »  -  «  for  0^e<l. 

If  the   candidate  list   is  empty,   stop. 
LB(e)   =  g(0)    for  o^esi. 

Select   a  candidate  problem  and  remove  it  from  the  candidate 
list.      Call   it    (CP^). 

Solve   (CP^  ).      If  it  la  infeasible,   obtain  the  appropriate 
dual  extreme  point    (ii*,   v*)    and  extreme  ray   (y*,    z*) . 
Otherwise  obtain  an  optimal  dual  solution   (u*,   v*) . 

I.  (CP^  )    infeasible. 

K 

(a)  y*d  ^0.   Go  to  Step  2. 

(b)  y*d  >  0.   Set  6*  -  [-y*(b-0*^)-z*e]  /  y*d. 

If  UB^(e;*)  <  LB(e)  for  all  6*20^1,  go  to  Step  2. 

II.  (CP^  )  feasible. 

If  UB^(e;*)    <.  LB(6)    for  all  059^1,    go   to  Step   2. 

If   the  optimal  primal  solution  of   (CP^  )    is   all  integer, 
use   it   to  update  LB(9). 

Use  a  heuristic  to  find  feasible  solutions  of   (CP^)   with 
right-hand-side   (b+Bd)    for  several  values   of  6.      Use   these 
feasible  solutions   to   update  LB(e). 

Separate   (CP'')    into  two  new  candidate  problems   (CP     )    and 

(CP^     )   by   choosing  p«  J*^  and  setting  F*^     -   f'^       -   f''   u   ip) 

x^     =0,   x^       -1.      Place   (CP'^  )    and   (CP^     )    In   the  candidate 
P  P 


list  and  return  to  Step  2. 

I 
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The  validity  of  the  generalized  bounding  test  insures   that   an 
optimal  solution   for  every   (P  )  ,   O^se^l,  will  be  found  by  the  search.      At 
worst,   an  optimal  solution  may  not  be  discovered  until  the  bottom  of  the 
branch-and-bound  tree  is   reached  (f'^-J).     This   guarantees   that  LB(e)   will 
coincide  with  g(e)   by  the  time  the  algorithm  is   finished.      It   remains  only 
to  show  how  the  optimal  solutions  are  identified. 

Let   {x   |keK}  be  the  set  of  incumbents  when  the   algorithm 
terminates.      Let   Be [0,1]    and  suppose  that   (P.)    is    feasible,   g(e)    >  -  «. 
From  the  construction  of  LB(6)  ,   (3.1)   -   (3.A)  ,  we  know  that   there  is  some 
keK     such   that 

g(e)    =  LB(6) 
=LB^(e) 


k 

I       r  X     >  -  "  . 
j=l     ^  ^ 


k  k  k 

Furthermore,   LB    (6)    >  -  <»     means   that   6^    26^   Q„    ,   or  equivalently   that 


n 

I       Ax 

j=l     ^ 


^  2  b+ed. 


Since  x  is  feasible  for  (P„)  and  its  return  is  equal  to  g(0),  it  follows 

o 

k  -  k 

that  X  is  optimal  for  (P„).  To  summarize,  if  keK  and  6e[0,l],  then  x  is 

0 

optimal  for  (P„)  if  and  only  if 

n      , 
i)    I       A  x^  <.   b+ed 
j=l  J  ^ 


and   ii)   I       r.x;  =  g(e)  . 
j=l  ^  ^ 


-12- 


At  Step  6,  in  contrast  to  the  prototype  algorithm,  x  is  not  fathomed 
when  the  optimal  primal  solution  of  (CP^  )  is  all  integer.   This  is  because 
X  may  have  other  descendants  which  are  optimal  for  6>0.   The  use  of  heuristics 
at  Step  7,  while  in  principle  optional,  is  an  Important  part  of  the  algorithm 
since  Integer  solutions  of  (CP^  )  can  only  yield  LB(e)  -  LB(0)  for  O^e^l. 
The  heuristics  are  needed  to  produce  stronger  values  of  LB(e)  for  6  >  0. 

As  with  the  prototype  algorithm,  the  above  procedure  will  admit 
considerable  variation  and  refinement.   If  the  dual  simplex  method  is 
used,  then  suboptimal  dual  solutions  can  be  used  to  perform  additional 
bounding  tests.   Cutting  planes  can  be  generated  for  any  candidate 
problem  to  give  stronger  upper  bound  functions.   Parametric  linear  programming 
can  be  used  to  generate  more  than  the  first  segment  of  UB^(9).   If  a 
candidate  problem  with  an  all  -integer  LP  solution  has  to  be  separated  at 
Step  8,  then  the  same  LP  solution  is  optimal  for  one  of  the  two  new 
candidates  and  does  not  have  to  be  recomputed.   Extensive  experimentation 
will  be  required  to  determine  the  most  effective  computational  tactics. 
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6.     Example 

In  this  section  the  algorithm  will  be  applied  to  a  simple  example. 
In  order  to  Illustrate  all  of  the  different  cases  that  can  arise,  the 
parameterization  will  be  done  over  a  relatively  large  interval.  The  test 
problem  Is 


max  lOx,  +  15x-  +  lOx,  +  5x, 
12      3     4 


subject  to 

2x,  +  3x,  +  5x,  +  Ix,  S  4  +  94 

12     3     4 

4x,  +  2x^  +  Ix^  +  Ix,  ^  4  +  64 
12     3     4 

X  e  {0,1}       Ujs4 

Thus  b=(4,4),  d=(4,4)  and  increasing  9  from  zero  to  one  amounts  to  doubling 
the  right-hand-side.  A  picture  of  the  optimal  return  function  fCb"*)  is 
given  in  Figure  3.   The  dashed  line  indicates  the  line  segment  of  interest: 
{b+6d  I  0S9S1}.   It  is  clear  from  this  picture  that  three  optimal  solutions 
must  be  found,  with  values  of  20,  25,  and  30.   These  solutions  are  (0,1,0,1), 
(1,1,0,0),  and  (1,1,0,1)  respectively.  The  g(e)  function,  shown  in  Figure  4,  is 


g(9) 


,  20  for  0  s  9  <  1/2 

25  for  1/2  ^  6  <  3/4 
I  30  for  3/4^0  2  1. 


The  optimal  LP  solution  of  (P  )  is  x=(l/2,l,0,0) ,  u-(5,0),  v-(0) 
with  value  UB  =20.   The  rounded  down  solution  has  value  15  and  is  feasible 
for  e>0;  the  rounded  up  solution  has  value  25  and  is  feasible  for  6^1/2. 
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Thls  provides  an  initial  lower  bound  function: 

LB(e)  =  rl5  for  0  ^  e  <  1/2 
.25  for  1/2  <  6  ^  1. 

The  complete  branch-and-bound  tree  is  displayed  in  Figure  5.   The  nodes 
will  be  discussed  in  the  order  in  which  they  were  created. 

Node  1.  LP  solution:  x=(0, 1,0,1),  u=(5,0),  v=(0)  ,  UB"'"=20.  UB"''(e;*)  =206+20. 
The  LP  solution  is  all  integer  and  is  feasible  for  9^0.  Therefore  the  lower 
bound  function  may  be  improved: 


LB(e) 


20  for  0  ^  6  <  1/2 
.25  for  1/2  S  e  ^  1. 


I 


The  bounding  test  for  node  1  is  shown  in  Figure  6.   Node  1  is  not  fathomed. 


Node  2.   LP  solution:  x=(l, 0,0,0),  u=(0,10),  v=(0) ,  UB^=10.   UB^(e;*)=40e+10. 
The  bounding  test,  shown  in  Figure  7,  is  not  successful.   Notice  that  if  we 

were  only  interested  in  solving  (P^)  we  would  be  finished.   Node  1  has  an 

2 
all  integer  solution  with  value  20  and  node  2  has  upper  bound  UB  =10<20-LB(0) . 


3 
Node  3.  LP  solution:  x=>(0,0,3/5,l)  ,  u-(2,0)  ,  v-(0,0,0,3),  UB  -11. 

3 

UB  (e;*)-8e+ll.   The  bounding  test,  shown  in  Figure  8,  is  successful  and 

node  3  is  fathomed. 


Node  4.   Same  as  node  1,  since  optimal  LP  solution  at  node  1  has  x^  =  1. 
Node  5.   Same  as  node  2,  since  optimal  LP  solution  at  node  2  has  x„  =  0. 
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Node  6.   LP  is  infeasible.   The  dual  extreme  point  is  u=(0,10),  v=(0)  and 

the  extreme  ray  is  y=(0,l),  z=(0) .   The  critical  value  of  6  is  (-y(b-6  )-ze)/yd-l/2. 

Thus  UB  (6)=  -  "  for  0^e<l/2  and  UB  (e;*)=40e+5  for  l/2^e<l.   The  bounding 

test  is  shovm  in  Figure  9. 


Node  7.   Same  as  nodes  1  and  4,  since  optimal  LP  solution  for  node  4  has  x  -0. 


Node  8.   LP  is  infeasible.   The  dual  extreme  point  is  u=(5,0),  v=(0)  and 

Q 

the  extreme  ray  is  y=(l,0),  z=(0).   The  critical  value  of  G  is  (-y(b-B  )-ze)/yd-l, 

Q  1, 

so  UB  (6)  =  -  °°  on  0<e^l  and  node  8  is  fathomed.  , 


i 


Node  9.   Same  as  nodes  2  and  5,  since  optimal  LP  solution  for  node  5  has  x  =0. 


Node  10.  LP  is  infeasible.  The  dual  extreme  point  is  u=(5,0),  v-(0)  and           ' 

the  extreme  ray  is  y"(l,0),  z=(0) .   The  critical  value  of  9  is  (-y(b-3  )-ze) /yd«='3/4. 

Thus  UB'''°(e)=  -  "  for  0^e<3/4  and  UB-'-^CG;*) -209+5  for  3/4^9^1.  Node  10  is 
therefore  fathomed.   See  Figure  10, 


Node  11.   LP  is  infeasible.   The  dual  extreme  point  is  u"(0,5),  v=(0)  and 

the  extreme  ray  is  y=(0,l),  z=(0).   The  critical  value  of  9  is  (-y(b-e  )-ze) /yd-1/2. 

Thus  UB''"''"(9)  =  -  <»  for  0^9<l/2  and  UB"'"''" (9;*)  =209+15  for  1/2^6:^1.   Node  11  is 

not  fathomed.   See  Figure  11. 


Node  12.   LP  is  infeasible.   The  dual  extreme  point  is  u-(5,0),  v-(0)  and 

12 
the  extreme  ray  is  y=(l,0) ,  z=(0) .   The  critical  value  of  G  is  (-y(b-3  )-ze)/yd-i 

12 
Therefore  UB   (9)  =  -  «  for  0^6^ 1  and  node  12  is  fathomed. 
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Nodes  13-18  are  all  at  the  bottom  level  of  the  search  tree.   The 
solution  for  node  18,  (1,1,0,1),  has  value  30  and  is  feasible  for  6^3/4. 
The  lower  bound  function  may  be  improved  by  redefining  LB(e)=30  for  3/4^9^1. 
LB(e)  now  coincides  with  g(e)  on  0<e^l.   The  algorithm  terminates  since  the 
candidate  list  is  empty. 


The  amount  of  extra  computation  required  to  solve  (PIP)  ,  as  compared 
to  (P  ) ,  depends  on  the  length  of  the  interval  of  parameterization.   When 
this  interval  is  small,  the  burden  imposed  by  parameterization  may  be 
slight  or  even  negligible.  When  it  is  large,  however,  as  illustrated  in 
this  example,  the  burden  can  be  quite  substantial. 


I 


"2' 


8 
7. 


4 
3 
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Figure  3.   The  optimal  return  function  f(b'). 
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Figure  4.   The  parametric  function  g(e) . 
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Flgure  5.  Branch-and-bound  tree  for  the  example. 
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30 

25 
20 

15 

10 

5 

0 


206+20,^-^^ 

LB(e) 

^^^^^ 

1/4 


1/2 


3/A 


Figure  6.  Bounding  test  for  node  1. 


30 
25 

20 
15 

10 

5 

0 


409+10    y 

LB  (9) 

X 

1/4 


1/2 


3/4 


Figure  7.   Bounding  test  for  node  2, 
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Figure  8,      Eoiindlng  test   for  node  3, 


I 


30 

400+5, 

9^1/2/^ 

y^            LB(e) 

20 

15, 

10. 
5 

1/4 


1/2 


3/4 


Figure  9.   Bounding  test  for  node  6. 
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Figure  10.  Bounding  test  for  node  10. 
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Figure  11.   Bounding  test  for  node  11. 
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7.   Computational  Results 

The  ideas  presented  above  were  tested  by  incorporating  them  into  a 
branch-and-bound  computer  code  [3].   The  results  for  three  test  problems  are 
presented  in  Table  1.   In  each  run  the  direction  vector  d  was  taken  as  some 
percentage  of  the  right-hand-side  b.  For  example,  if  d-5%b,  then  (PIP)  has 
right-hand-sides  b+e(.05)b  for  Osesl.  The  column  headed  "solutions"  gives 
the  number  of  optimal  solutions  found,  or  equivalently  the  number  of  steps 
of  the  g(e)  function.   "Heuristic"  is  the  number  of  (evenly  spaced)  6  values 
for  which  the  heuristic  is  applied  at  Step  7.   The  problems  are  of  the  capital 
budgeting  type  and  the  heuristic  employed  is  that  of  Toyoda  [11].   "Pivots" 
is  the  total  number  of  linear  programming  pivots  and  "time"  is  the  total 
solution  time  in  seconds  on  an  IBM  370/168. 

These  results  illustrate  quite  clearly  how  the  computational  burden 
increases  as  the  interval  of  parameterization  is  lengthened.   In  order  to 
facilitate  comparison  with  our  results  by  other  researchers  we  have  included 
the  data  for  the  5x30  problem  as  Table  2  and  the  corresponding  g(9)  function 
for  a  10%  increase  in  b  as  Table  3. 
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Table  1.     Computational  Results   for  three  test  problems. 


m 

n 

d 

solutions 

heuristic 

pivots 

time 

5 

15 

0 

1 

1 

39 

.239 

.05b 

4 

10 

62 

.541 

»10b 

5 

10 

91 

.815 

.15b 

7 

10 

124 

1.044 

.2Gb 

8 

10 

130 

1.170 

.25b 

10 

10 

171 

1.534 

.50b 

16 

20 

315 

3.162 

5 

30 

0 

1 

1 

153 

1.605 

.05b 

10 

10 

479 

7.280 

.10b 

26 

10 

1958 

24.040 

10 

28 

0 

X 

1 

67 

1.073 

.05b 

16 

10 

181 

4.192 

.10b 

29 

20 

637 

11.970 

.15b 

42 

20 

12  &6 

29.204 
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Table  2.      The  5x30   test  problem. 


^iJ 


zj 


*3j 


*4j 


'5j 


b  = 


188 

91 

20 

86 

164 

936 

92 

179 

99 

97 

98 

695 

6 

146 

95 

42 

2 

390 

80 

155 

95 

90 

165 

1152 

91 

102 

84 

101 

140 

980 

44 

112 

136 

3 

106 

1000 

108 

126 

166 

101 

88 

815 

166 

21 

13 

34 

68 

109 

171 

39 

20 

25 

84 

807 

64 

67 

124 

72 

131 

156 

97 

29 

42 

96 

55 

548 

35 

55 

58 

36 

11 

335 

51 

72 

43 

3 

17 

316 

98 

17 

43 

88 

4 

528 

36 

0 

44 

97 

47 

36 

70 

42 

2 

77 

45 

573 

27 

15 

88 

50 

11 

38 

94 

64 

55 

14 

77 

3 

68 

53 

68 

77 

36 

800 

13 

30 

22 

88 

49 

392 

13.2 

2,8 

6.8 

11.3 

2.9 

92 

15.1 

15.0 

8.3 

13.8 

11.7 

4 

3.3 

2.6 

8.9 

4.5 

19.2 

29 

7.4 

3.5 

3.1 

17.1 

18.1 

81 

7.0 

17.0 

16.5 

11.8 

3.8 

2 

1.2 

3.5 

2.2 

17.1 

18.0 

40 

7.0 

5.1 

9.7 

19.1 

8.8 

17 

17.0 

16.2 

4.7 

5.0 

3.9 

16 

13.8 

13.2 

1.8 

10.2 

16.9 

30 

9.4 

13.9 

11.0 

3.6 

13.8 

118 

800 

800 

700 

700 

800 
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Table  3.      The  g(e)    function   for  a  10%  Increase  in  b;   5x30  probl 


em. 


e 

g(e) 

e 

g(e) 

e 

g(9) 

0.0 

7515 

.A4166 

7806 

.74142 

7947 

.02750 

7578 

.52499 

7822 

.74714 

7957 

.09666 

7607 

.59000 

7839 

.78499 

7987 

.15428 

7612 

.64250 

7849 

.81285 

8009 

.16875 

7633 

.64285 

7850 

.84428 

8049 

.17125 

7696 

.65571 

7891 

.89124 

8060 

.20374 

7725 

.68500 

7913 

.93832 

8079 

.33625 

7726 

.70833 

7931 

.96666 

8112 

.36666 

7777 

.71750 

7942 
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